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■ Abstract 

We compute, on the (A$ )i+i model on the lattice, the soliton 
mass by means of two very different numerical methods. First, we 
make use of a "creation operator" formalism, measuring the decay 
' of a certain correlation function. On the other hand we measure the 

shift of the vacuum energy between the symmetric and the antiperiodic 
systems. The obtained results are fully compatible. 

We compute the continuum limit of the mass from the perturbative 
Renormalization Group equations. Special attention is paid to ensure 
that we are working on the scaling region, where physical quantities 
remain unchanged along any Renormalization Group Trajectory. We 
compare the continuum value of the soliton mass with its perturbative 
value up to one loop calculation. Both quantities show a quite satis- 
factory agreement. The first is slightly bigger than the perturbative 
one; this may be due to the contributions of higher order corrections. 
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1 Introduction 



Standard perturbation is known to be a useful tool for the formulation of 
Quantum Field Theory starting from Classical Field Theory. It has, however, 
serious handicaps such as the fact that non-perturbative effects are not taken 
into account. An alternative possibility is the quantisation of non-trivial, 
non-perturbative solutions to the classical equations, such as solitons. 

The study of such topologically non trivial vacua in Field Theories pre- 
sents several problems when the model is formulated in the continuum or in 
the lattice. In the continuum it is very difficult to extract non perturbative 
quantities, like mass, and in the lattice, where this is possible with Monte 
Carlo simulations, other problems are present. 

First, on the lattice the analysis of this kind of configurations is made 
difficult by the trivial topology of the lattice, and because the concept of 
continuity is lost Second, it is of fundamental importance to define how 
to measure quantities on these topologically non trivial vacua 0. We can 
consider how to compute the mass, for instance. As is known, solitons are 
characterized by a topological charge, related to their behaviour when spatial 
coordinates tend to infinity ( Q ~ $(x oo) — $(x — oo) ). This charge 
is conserved with time. When we quantise a soliton, we obtain a quantum 
soliton-particle and a series of excitations of this particle, a so called soliton 
sector. Topological charge becomes then a quantum number characterizing 
the sector. Its conservation prevents the soliton from falling to the vacuum, 
ensuring its stability. The standard way of calculating the mass is considering 
an operator with non vanishing projection on this sector, then computing the 
connected correlation to large distance, and finally extracting the mass from 
the coefficient of the exponential decay. 

In the general case, for topologically non trivial sectors, the definition of 
such an operator is very ambiguous. It is possible to define many operators 
on the lattice sharing the same continuum limit, although their behaviour 
far from this limit differs from each other. On the lattice, the region where 
we can obtain results within reasonable computation times is generally far 
from the continuum limit -where very big sizes would be necessary-, and all 
those continuum-equivalent operators give us different results, 

In four dimensional theories computer limitations have made this point 
particularly difficult. 

Fortunately, some interesting facts can be studied quite satisfactorily in 
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less than four dimensions. We consider on this paper the (A0'^)i+i model, 
where solitons are also present. 

On a finite lattice the boundary conditions fix up the topological sector. 
Periodic conditions fix the trivial vacua, for instance. Antiperiodic conditions 
fix vacua with non trivial topology (if the symmetry is broken). Only free 
boundary conditions allow us to have different topological sectors; however 
the finiteness of the system allows us to travel between vacua, and we finish 
always in the trivial sector, the energy of which is lower. 

In this model it is possible to carry out the computation of the soliton 
mass by using two different, related, methods. 

First, we have made use of the operator defined by Kadanoff et al. [|]; 
in Spin Systems, its effect can be seen as the introduction of a twist: a 
topological excitation induced by a specific dislocation of the lattice. It has 
a topological charge different from zero. Consequently, we expect a non-zero 
projection onto the sohton sector. 

On the other hand, we can consider the system with antiperiodic spatial 
conditions for the scalar field. This system can be considered as the periodic 
one after the introduction of a twist along the whole lattice time. The differ- 
ence of the energies of the periodic and antiperiodic systems, which is a local, 
easy to measure quantity, provides us with another method to compute the 
soliton mass. 

We always keep in mind that a theory in a lattice acquires physical mean- 
ing only when we make its spacing tend to zero. In order to get to the con- 
tinuum limit 0] we use the Renormalization Group (RG) equations, which 
are known in this model. We must consider the limit of zero lattice spacing. 
A change in this spacing, and a change in the coupling constants in such a 
way that the physical observables remain unchanged, can be carried out by 
using the RG equations. 

Iterating RG transformations, we obtain a series of points in the pa- 
rameter space, - Renormalization Group Trajectory (RGT)-, which can be 
characterized by a parameter /. Different points of a trajectory, correspond- 
ing to different values of the parameters, are obtained after integrating over 
successive energy scales. Thus, when we move on any RGT, the Physics 
remains the same. In this way, for example, as we evolve on a RGT, we 
see different values of the correlation length on lattice units, ^o{^), but this 
correlation in physical units ^ = a{l) ■ ^o(0 remains constant. 
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1.1 model 

We study the model in d=2 dimensions, the euchdean lagrangian density 
of which is given by 

C-... = \{d,^f + "^^' + ^^' , (1) 

where $ is dimensionless, [A] = [r]= 

In order to adapt our lagrangian to the lattice, we proceed as usual: 

n=0 ^ 



where L is the lattice extension. We conclude 

2^ " ' 4! 



5euc = - E '^n'^'n.M + Eil^ + + ^<^n} , (3) 



n,fi 



We introduce the following notation: the adimensional parameters defined on 
the lattice are subscripted; thus we use Aq, tq (respectively equal to Aa^, ra^). 

Making the spatial coordinates discrete implies imposing a momentum 
cut-off A = After scaling the momenta g — p = we can express in 
momentum space H: 



Z Jp 4! Jpi Jp2 Jp3 

(4) 

with / = Jq ^ In d=2 we have two fixed points 0|: 

i) The gaussian point, S'gauss = /p ^2(p)'^(p)'^(~P)' with U2{p) ~ p^. That 
is to say, taking just the kinetic part of the lagrangian. 

ii) A non-trivial point, which is built adding to the lagrangian the term 

S= uliii^(i2,(i3,-qi-q2-q3)^{qi)^iq2)^{q3)^{-qi-q2-q3) , 

with ul{qi...qi) ~ (g? + ... + 



The A$^ model in less than four dimensions is superrenormalizable. The 
only divergent graph is that of one vertex with two external legs and a loop. 
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We can get rid of this divergence simply by renormalizing the mass, and 
therefore it is not necessary to renormahze A; we can keep it fixed all the 
time as we make a go to (equivalently, A — > oo). Since A = Aq a~^, it 
implies Aq as a^. In our lattice, consequently, in the continuum limit 
U4 = Xq = 0'. we are considering the gaussian fixed point. Our RGT will 
evolve in its attraction domain. 

We follow the Renormalization Group scheme; in order to permit a conti- 
nuous evolution in the parameter space, we allow integrations of the variable 
p between ^ and 1. Linearising the resulting equations near the critical point, 
we obtain |^ 

Xo = s^"'X,; fo = s2{ro + ^logs}, (6) 

where Aq and fg are the transforms of Aq and r^. 

These expressions have a limited region of validity: for big values of Aq, tq, 
the linear approximation is not valid. On the other hand, for small values of 
the parameters, we are very near the gaussian point, the correlation length 
becomes very large, and if it is of the order of the lattice size, finite size effects 
mask our results. We refer to the intermediate region where the continuum 
is reproduced as the scaling region. 

We remark that the fact of staying in the basin of attraction of a Gaussian 
fixed point does not prevent at all the possibility of spontaneous symmetry 
breaking. Given one point (AcTq) in the parameter space, the Renormaliza- 
tion Group Trajectory starting from it cannot cross the transition line be- 
tween the ($^) = and ($^) 7^ phases; it remains in the phase to which 
the initial point belongs. Thus, if we start in the symmetry broken phase, 
the continuum limit of our theory presents symmetry breaking. 

2 Computation of the soliton and fundamen- 
tal boson masses 

2.1 Fundamental boson 

In order to calculate the mass of the fundamental boson rup we use the 
connected correlation function between the Higgs fields, 0)$(x, t)). In 
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order to avoid contributions from states with non-zero momenta, we integrate 
on X and consider 

c^it) = imm), (7) 

where = / dx^{x,t). 

For large t, and if the correlation length is different from zero, C^{t) 
behaves, on an infinite lattice, as C^{t) ~ exp(— mt). 

We consider periodic boundary conditions in our lattice, L being its extent. 
Consequently, the point n is equivalent to n + L. Given two points at a 
distance t, there are two possible paths connecting them: one of length t 
and the other, resulting from the boundary conditions, the length of which 
is L — t. Thus the mass is given by 

(0(t)0(O)) ~ e-"* + e-"^(^-*) (8) 

so that 

C^{n+ 1) _ cosh[mp{n + 1 - |)] 
C^{n) ~ cosh[mp(n + -f )] ' 

where a is the time spacing of the lattice. We can solve (P) and obtain a series 
of values of mp{n) depending on n. For small n they have contributions from 
large mass states, and for large n the signal is small; there is an intermediate 
region of n where mp{n) is nearly constant. We take it as the actual value of 
the mass. 



2.2 Soliton mass 

Kadanoff introduces the correlation function between two points in the 
dual space i?2, {f^Rif^R2) iii the following way: we start from a lagrangian 
S = J2n,/i Jn,fi^n^n+fj, + J2n 0{^n)', wc draw a path in the dual space connect- 
ing the dual points -Ri and R2, and change the sign of the coupling constants 
J's placed on the links crossed by our path. We have 

ifiR.fiR^) = ^^exp{-^*J„,^$„$„,^ + ^'J„,,^$„$„,^}, (10) 

where J2[<s>] runs over all the configurations of the field, Yln,^l takes into ac- 
count the links with their signs changed and J2'n fi refers to the rest of the 
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links. Equivalently, 



[<I>] "lA* 'T'lM 

(11) 

where S is the original action. 

We can alternatively express the correlation function in terms of a "twisted 
action" St = 5* — 2 ■ J2n,fi Jn,fj.^n^n,iJ., with it corresponding partition function 
Zt = J^m 6xp{— Sj}, that is to say, 

In our case, with a lagrangian given by (0), we can express it in a similar way, 
depending on link variables J„,^, after making a change of variable $ — ^ 
Now, J„^^ = constant = J > 0. This causes the appearance of a twist: the 
fields placed in the points n, n + fi connected by a link where J„,^ has changed 
to - Jn,^_i tend to change their signs: we have given rise to a topological 
excitation, with a non-zero topological charge $(x = oo) — $(x = —oo). 
We now define 

C^(^) = (Mr^,t)M« + r)), (13) 

where our path in the dual lattice will be the minimum length path connect- 
ing them, i.e. straight vertical lines. 

The topological excitation (with non-zero projection on the soliton sec- 
tor) appears at the time t, and annihilates at t + r. Thus, we expect an 
exponential behaviour, similar to that of C^{t). In this case, one of the two 
paths connecting the points i?i and R2 has a much bigger contribution to 
C^: that which crosses the dislocation. In fact, we have observed a clear 
exponential decay, (/i(n:,t)/U(ri,t+T)) ~ exp{— Cr}, and therefore we obtain the 
soliton mass as 

m.„.... = -log^-^. (14) 

When we study (e"^^"-'' we must consider the risks of our method: 

we study a strongly non-local quantity, which is seriously affected by the fi- 
nite size of our lattice. Besides, the use of an exponential function implies a 
magnification of errors. 
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In principle, we do not know to what point these effects will spoil our 
results. In order to control these risks, we look for an alternative way of 
calculating the soliton mass from local non-exponential variables. Following 
Groeneveld et al. we introduce a local parameter which accounts 

for the energy response to the appearance of the twist. 

First, it will be useful to change the variables the action depends on. 
We note that, making the following change of variable, $ — > ^ = —y=- with 



^0 



/3o = ^, we obtain 

^(ro, Ao) = Z{ro,Po) = 

_v f+oo J- 1 

where V = L'^ is the volume of the system. The (3q derivative of Z{rQ, Aq) 
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dZ{r„,l3„) V Z{r„,i3„) 1 

where S'(ro, Aq) is the action resulting from the integration of our lagrangian 
(i- 

If we impose antiperiodic boundary conditions in the spatial direction, 
we introduce a twist the length of which is the temporal dimension of the 
lattice, T. We can define the "twisted" partition function corresponding to 
this twist, Ztir^X) = Zt{r,P). Keeping in mind (|I2|) we can now calculate 
the soliton mass as 

m^oi = --^ log(Ai(n, T)i2{n, 0)) = log - 



where jSc is the value for which our trajectory cuts the transition line between 
the ($^) = and ($^) ^ phases, and 

m) = -^7^1og§r^ = ^kSt{ro,Xo))t - (5(ro,Ao)) , (18) 
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where () and {)t stand for expectation values with Z (periodic boundary 
conditions) and Zt (twisted or antiperiodic) respectively. We remark that 
in the Zt system ($) = in both the symmetric and broken phases. The 
integration in ( p!7D implies defining a trajectory in the parameter space with 
ro fixed, and /3 starting from (3^- For higher values of (3, we have no symmetry 
breaking, and the soliton mass vanishes. 

We will check the masses obtained with the exponential function by com- 
paring them with those resulting from using Q. 

3 Details of the simulation 

We have made use of a specially-designed transputer based parallel machine, 
RTN, including 64 T-805 processors distributed in 8 boards with 8 each. 
As an individual board calculates one point in the parameter space, we get 
eight absolutely independent groups of measurements for every (AcTq). The 
error for every magnitude has been calculated averaging its 8 independent 
predictions. We have used an adaptative MC process so as to keep the rate 
of acceptance between 40 % and 60 %. 

We have simulated different lattice sizes (16^ 24^, 48^ ), with (1000, 3000, 
7000) iterations of thermalization and (2000, 22500, 30000) measurements. 
Within each transputer, we have taken (20, 10, 5) decorrelation MC itera- 
tions between two consecutive measurements. We have observed no relevant 
finite-size effects in the local quantities. However, big sizes are needed when 
computing correlations, especially those defined by Kadanoff's operator, as a 
consequence of its strong non-locality and its exponential form. As we have 
mentioned earlier, we can not rely on the small-length correlations because 
of the contribution of large mass states; on the other hand, long distance 
correlations are seriously affected by the finiteness of our lattice. Thus, in 
order to obtain a precise value for the mass from Kadanoff's operator we 
have needed bigger and bigger lattices. 

4 Phase diagram and the scaling region 

Our model exhibits two phases. Classically, for positive values of tq the 
minimum energy configuration is $ = 0; for negative tq, a spontaneous 
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symmetry breaking occurs, and the new minima are $ = iy -j^. 

When we consider the contributions of all the configurations, each weighted 
with exp{—Seuc}, for small negative values of tq both minima are very close 
to each other and are not deep enough to stop fiuctuations from restoring 
the symmetry. More negative values of tq are necessary to ensure that we 
are in the broken phase. Therefore, in the semiplane with negative tq there 
is a transition line separating both phases. 

In order to determine the transition line, we choose several values of Aq. 
For each of them, we decrease ro until (<l'^) becomes different from zero; in the 
limit of an infinite volume, its value passes from zero to a finite non-zero value 
when crossing the transition line; in a finite- volume system in the symmetric 
phase, (<l>^) ~ and what we see is a sharp rise of ($^) (technically, in 

a finite lattice ($) is not a good order parameter because tunnelling between 
states with positive and negative values of the field cause it to be equal 
to zero all over the parameter space). Another useful quantity as an order 
parameter is the soliton mass. When computing {exp{—2JJ2^n^n+fi}), if 
we are in the ($^) = phase, the values of fiuctuate around zero, and their 
sum over the path vanishes; the expected values appearing in become 
independent of the length of the path and equal to 1 and rrisoi is zero. On the 
other hand, we expect non-zero nisoi for the symmetry broken phase, where 
the J2 ^n^n+ii is different form 0. We determined the transition line using 
both parameters ( ($^) and mgoi)- It is shown in figure 1. 

In the region where ^ is large enough, we can compare them with mean 

field predictions. Theoretically, |$| = \/^^- In this region the fiuctua- 
tions are small, and thus we have {exp{—2JJ2L ^n^n+ij,}) ~ exp{— 2 JL|<I>^|}, 
where L is the length of the summation path, and we can consider ($(0)$(r)) ~ 
|$^| and therefore rrisoi = 2J|$^|. Our results agree with these predictions. 

Now we pass to determine the scaling region. We keep ourselves in the 
(<l>^) 7^ phase, where rup and rusoi are different from zero. The reason for 
this is that, proceeding in this way, the continuum limit of our theory will 
correspond to the symmetry broken phase of the continuum problem, which 
is the one we are interested in. 

Thus, our next step is finding the region in the parameter space where 
equations (|^) are valid. 

Along a RGT we expect to find constant values for the physical mean- 
ingful variables, such as the correlation length ^, the physical masses M... In 
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the lattice we work with adimensional quantities depending on the point 
of the trajectory {^o{ro, XQ),m{rQ, Xq)...), related to the physical ones by 
M = = a ■ ^o--- Consequently, although our lattice-defined quantities 
vary, the ratios mp/nisoi , fn'^p/Xo , fn'ioiiton/ ^o, which are equal to the phys- 
ical expressions Mp/Mgoi , M'^/X , M^^JX, remain constant along these 
trajectories. 

We start from different points in the parameter space (rg, Aq) near to the 
critical point (—0.2 < ro < 0, Aq = fixed = 0.1). Iterating Renormalization 
Group transformations (eq. (P), where we choose s=1.08), we get further and 
further away from the origin (and thus from the continuum limit), drawing a 
series of trajectories. Along each one of them, we calculate the previously de- 
fined ratios in the different points obtained by the transformations. For every 
trajectory, we find a segment where these quantities remain approximately 
constant; the union of all the segments gives us the scaling region. 

Initially, we follow curves near the transition line separating the ($^) = 
and ($^) 7^ phases. As we move away from it, we find that the length of the 
segment reduces. This is clear, because we need a large lattice correlation 
length in order to reproduce the continuum limit, and the region close to 
the line transition is appropriate to that. Far from this line the correlation 
length is small, and the discretization is important. 

Finally, we choose a curve near that line, with the initial values 
(ro= -0.105, Ao= 0.25) (see figure 1). At this point we can illustrate our 
comments about the difficulties derived from the use of Kadanoff 's operator. 
In figure 2 we represent the correlations of the fields $ and /i for some points 
of this trajectory. As expected, when we get near the gaussian point, the 
correlation length increases, the mass is lower and the correlation decreases 
more and more slowly. For small enough values of the parameters, the cor- 
relation function {iiNl^N+n), for distances of the order of the length of the 
lattice, is not compatible with zero, (n is the distance in units of the lattice 
spacing). Thus, we must be very careful when we calculate masses in this 
region. 

There is another reason that makes it desirable to work with big lattices. 
Our method for calculating the masses consists basically on finding a certain 
correlation function, and fitting it to an exponential, or to an hyperbolic 
cosine. We expect this fit to be reasonably good for a set of intermediate 
values of n. When n approaches the length of the lattice -in our case, half this 
length, because of the periodic boundary conditions-, the fit is not possible 
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any longer. The bigger our lattice is, the longer this well-fitting segment 
becomes, and we have more points to fit our theoretically predited behaviour, 
and so calculate the mass with higher precision. This is clearly shown in 
figure 3: in the plot at left we draw the logarithm of the correlation function 
{^NfJ'N+n) for 24^ and 48^ lattices, in a region far from the gaussian point. 
In the small lattice, when n ~ 9, the fitting to a straight line is no longer 
possible, while in the big one we can still include some more points and get 
a good fit to a straight line. In the small picture at right, the parameters 
are Aq = 0.25, vq = —0.105; we have seen in figure 2 that, for these values, 
the correlation length is comparable with the lattice length, and we expect 
serious corrections. In fact, for the small lattice, the agreement region is 
smaller. 

In figure 2 we see that the function {^n^n+u) is much smaller than the 
H correlation, and we expect that the values obtained for the boson mass are 
not so strongly affected by the size of the lattice. Our results confirm this 
prediction. 

Now we can estimate the scaling region. From figure 4, we see that it 
begins at Aq ~ 2, the value from which mp/msoi can be considered as a 
constant. The upper boundary of this region can be more clearly inferred 
from figure 5. We expect m/\/Xo to be constant or, equivalently, a linear 
behaviour of m with a/Aq, m tending to zero as \/\) does. Thus, for the 
RGT starting from the initial values Aq = 0.25, Tq = —1.05, the scaling 
region corresponds to the interval 1.5 < Aq < 14 

5 Results 

5.1 Results from the operators 

First of all, we want to check that what we call sohton mass, calculated using 
Kadanoff 's operator, really behaves as a mass. We will compare its evolution 
under the Rcnormalization Group equations to that of the fundamental bo- 
son mass. On the other hand, we compare its value to previous theoretical 
predictions. 

A quantity Qphys with dimension [©] is related to its equivalent in the 
lattice Qiatt by Qphys ~ fl^®'©iatt- As we are working in less than four di- 
mensions, we can keep A fixed, and Aq ~ A • tends to as when we 
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approach the continuum hmit (gaussian fixed point) along a RGT. We also 
expect finite values for Mp,Msoi, so our value of the masses in the lattice 
m M ■ a tends to as a. 



From that we deduce, as we approach the continuous limit, rusoi ~ \/Ao. 

Figure 5 shows that, within the limit of the scaling region, that is the be- 
haviour of nip and uisoi- 

In a system with only one relevant direction, the renormalized trajectory 
coincides with that direction. In our case we have a twice-unstable point, and 
there is a continuous family of renormalized trajectories leaving it, each of 
them corresponding to a different continuum theory. In the previous section 
we have chosen one of those trajectories; once we give an arbitrary value of A, 
Mp and M,oi can be calculated as M = -^^A. In our RGT, Mp = 0.453a/A 

and Msoi = 0.356VA. 

Next we study the evolution of rrisoi when, starting from a point on the 
transition line, we move further and further into the symmetry broken phase. 

First, let us summarise some qualitative basic ideas. When we quantise 
the classical absolute minimum, we obtain the vacuum of the quantum the- 
ory; quantisation around the soliton gives us the sohton sector. If this local 
minimum is broad (which, in our case, corresponds to a point near the tran- 
sition line between ($^) 7^ and ($^) = 0), a great number of configurations 
different from the classical solution will contribute to the value of any ob- 
servable. But as we move away from this zone, the potential well gets deeper 
and deeper, and a moment comes when we have contributions only from the 
minimum and configurations very close to it; we are recovering the classical 
solution. 

In order to explore these ideas, we trace a trajectory in (AcTq) space 
fixing ro= -2.2 and letting Aq move from the vicinity of the transition line 
(Ao — 12) deeper and deeper into the ($^) 7^ phase. This path cuts a 
different RGT in each of its points (with different physical masses for a value 



To calculate the classical continuous limit for the soliton mass, we make 
use of the fact that, in a continuous euclidean space, a soliton of the A^"^ 
classical theory in 1-1-1 dimensions propagating with a velocity v is given by 




of A). 




(19) 
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The energy density is e(x), the expression of which coincides with the 

euchdean lagrangian. The classical absolute minimum is $ = At a 

classical level, we can define the soliton mass as 



3 

r| 2 



Maass = E- Emin = 4^2^ , (20) 

Emin being the energy of the absolute minimum, and E that of the soliton 
solution. 

When we quantise the soliton, in the weak coupling approximation (— << 
1), the mass obtained is, up to an order 0(^A/|r|) [0 

1/3" 3 

Mquantum = M^ass + ^V^hx ^) • (^l) 

O V ^ TT 

We have taken all along h equal to 1. 

Comparing our results for rrisoi to those predicted by expression (^) (see, 
e.g., figure 7), we see a clear linear behaviour with jSo = I/Aq. That behaviour 
is intermediate between the 0(0) and 0{h) theoretical predictions, closer to 
this second one. That displacement with respect the 0{h) prediction may be 
attributed to the contribution of higher orders in ^Ao/|ro|. 



5.2 Results from twisted system 

The use of twists is known to be, in a computer simulated theory, a good 
help for studying the phase structure. Our main motivation for introducing 
it is to calculate rusoi in an alternative way to the use of Kadanoff 's operator, 
so avoiding its risks, already mentioned in section 2. 

We impose antiperiodic boundary conditions in the spatial coordinate. In 
this way, we introduce a twist in the lattice that lasts from t=0 to t=T. We 
evaluate the expected value of the action under these conditions, {St)t, fol- 
lowing the notation introduced in section 2. In the ($^) = phase, because 
of its ±$ symmetry, changing the signs of some J's does not cause the ex- 
pected value of the action to change, and (5*) = {St)t- That means rrisoi = 0, 
or we can also see the vacuum as a "soliton condensate". However, in the 
($^) 7^ phase things are not so any longer. By inducing the twist, we favour 
the appearance of a soliton propagating through time. 55" = {St)t — {S) is the 
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energetic response of the lattice to the introduction of the twist; it must be 
related to the euclidean energy of the soliton, as we will see. 

As we did in the previous section, we follow the path fixing the value of 
tq = -2.2. 6S must increase its value from near the transition line to the 
classically predicted one. Again, we plot 6S against ItqI^/Aq. In figure 6, 
we can see that it grows up steeply and soon stabilizes at the classical value 



(ID: \6S\^/Xo = 32|ro|7Ag. Instead of Aq we have drawn (3o = 1/A, 



0- 



Intuitively, there is a relationship between 6S and the soliton mass, which 
is its minimum energy level. In the classical limit, 6S/T is equal to rrisoi- In 
figure 6, we see that, as we increase Pq, both values tend to coincide. In 
general, the expression relating both quantities can be obtained if we keep 
in mind (^) and (|TB|) 

From fig. 6, we see that 6S/T suffers a sharp increase around /?o = 0.0826, 
(Ao = 12.1) indicating that we have crossed the transition line. We deduce 
f3c = 0.0804 ±0.0022 Ac = 12.45 ±0.35. In order to avoid errors coming from 
the estimate of jSc and 6S in the vicinity of that line, we take, instead of (^) 

msoiW) = m,,KA) ± ^ , (23) 

where we have taken /3j = 0.0826, the first value of (3 for which SS is clearly 
non-zero. We fit our values of nisoi to a straight line, and take mgoiiPi) to be 
the height of the line at f3i. 

We compare the results obtained applying (0) with those from Kada- 
noff's operator (see figure 7). Both values coincide with a precision up to 3 
% in the least favourable point. 



6 Conclusions 

We have studied topological excitations in the (A0'')i+i model on the lattice 
using a "disorder parameter", from the decay of which we can compute the 
soliton mass. The results obtained have a well defined continuum limit, which 
we have computed with the renormalization group equations. We have paid 
special attention to make sure that we are working in the scaling region, 
where physical quantities are unchanged along the RGT. 
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We have also computed the sohton mass from the difference of the vacuum 
energy between the twisted an untwisted systems, where quantities are local, 
and we found this result agrees with te previous one. 

Now we would like to compare both methods: the first has the disadvan- 
tages that it implies the calculus of operators which are strongly non-local 
and exponential, and thus the method is very sensitive to finite-size effects 
and systematic errors. A lot of statistics is necessary to obtain results within 
a reasonable margin of error. On the other hand, the method using twisted 
systems decreases considerably the computation time required. The statistic 
errors are small, and thus we conclude that imposing twisted conditions is a 
very satisfactory alternative in order to calculate the soliton mass. However, 
twisted conditions modify the vacuum of the theory, and make the calcu- 
lation of other masses to which to compare the results (such as, e.g, the 
fundamental boson mass) impossible. 

We have compared our non-perturbative result for the mass of the topo- 
logical excitations in the continuum limit with the theoretical perturbative 
result up to first order. Our results show a systematic lineal raise, which 
may be due to higher order corrections not considered in the perturbative 
prediction. 

The inclusion of fermions with a Yukawa coupling to the scalar fields is a 
very interesting future work. In this case we have a three parameter space and 
a very rich model. The problem is simplified by the fact that bosonization 
is possible, and so Monte Carlo simulation is simpler than when fermions 
are considered directly. Therefore, it must be also possible to compute the 
soliton mass and the continuum limit. 
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Figure captions 

Fig. 1 The Renormalization Group Trajectory followed starting from the ini- 
tial values To = —0.105, Aq = 0.25 is shown. The transition line between 
the ($^) = and ($^) ^ phases is drawn. 

Fig. 2 Decay of the correlation functions (/x//) and in a 48^ lattice. While 
the latter decays quickly to 0, {/ifi) is non-compatible with for large 
distances if we are near the gaussian point. Thus, we expect the value 
of the mass we calculate from it to improve as we increase the lattice 
length. 

Fig. 3 Logarithm of the correlation function (/i/i) for three different points 
in the parameter space: for the small window at right, Aq = 0.25, 
To = —0.105. The two other points are Aq = 2.934, tq = —0.945 and 
Ao = 5.431, ro = —1.616. Both are inside the scahng region. The latter 
is the one with the largest slope (the biggest mass) . 

Fig. 4 Evolution of the ratio nip/msoi along the RGT previously drawn in 
figure 1. From the results in the 48^ lattice, we see that the ratio 
becomes constant from Aq ~ 1.5, a fact which is not apparent in the 
24^ lattice. This gives us a lower limit for the section of our RGT inside 
the scaling region. 

Fig. 5 rrip and rrisoi versus ^/Aq along the RGT, for 24^ and 48^ lattices. We 
determine the scaling region by selecting the points which give a rea- 
sonable fit to a function y = a ■ x. In this way, we find the upper limit 
of this region, given by Aq ~ 14. For nip only the results from the 
48^ lattice are drawn because they coincide with those from the small 
lattice. 

Fig. 6 In a 48 X 48 lattice, along the ro = —2.2 vertical path, the evolution 
of the quantities {SS/T){Xl/\ro\^) and ?w^„^(Ao/|ro|^) versus Po = 
is shown. The classical value for both ratios is 32. The results are 
compared to those for the perturbative calculation up to an order 
0(Ao/|ro|). While all these results coincide in the limit (Ao/|ro|) — > 0, 
{6S/T) soon stabihzes at the classical value, while nisoi keeps closer to 
the first-order calculation. 
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7 For the 48 x 48 lattice, along the vertical path with tq = —2.2, the 
results for rrisoi obtained from the Kadanoff's operator and the mass 
given after the integration of fi(/3o) are compared. We compare them 
with the classical and the theoretical perturbative value up to first order 
in 0(Ao/|ro|), The x-co ordinate is Po — j-- 
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